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I.  Introduction 

This  program  Is  an  extension  of  the  work  performed  under  contract 
F30602-74-C-0082  and  reported  In  the  Final  Technical  Report  released  In 
January  1977  (RADC-TR-77-18).  A major  contribution  of  that  work  was  the 
exposition  of  a complete  theory  for  the  statistical  effects  of  atmospheric 
turbulence  on  coherent  radiation  reflected  from  a diffuse  target,  in- 
cluding the  effects  of  both  phase  and  amplitude  perturbations.  In  this 
report,  we  discuss  an  extension  of  the  previous  work  utilizing  the  complete 
theory  developed  there  and  including  a thorough  experimental  investigation. 

The  objective  of  the  present  work  was  to  extend  the  previous  results, 
both  experimentally  and  theoretically,  in  two  directions:  one  to  evaluate 
the  effect  of  target  glint,  and  the  other  to  consider  the  influence  of 
source  coherence.  Both  of  these  directions  have  led  to  significant  results 
that  will  impact  system  considerations. 

• With  a single  glint  on  an  otherwise  diffuse  target  we  have 
found  that  the  normalized  variance  of  the  scattered 
irradlance  is  enhanced  over  that  from  a diffuse  target 
alone.  This  increased  normalized  variance  (on  tlje  order 
of  3)  corresponds  to  a substantial  increase  in  the 
fluctuations  of  the  received  irradlance. 

• With  a small  but  finite  source  bandwidth,  we  have  found 
that  the  normalized  variance  of  the  irradlance  scattered 
from  a diffuse  target  is  unexpectedly  decreased  from  that 
when  a "coherent"  source  is  used.  This  corresponds  to 

a substantial  decrease  in  the  fluctuations  of  the  received 
Irradlance. 
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In  addressing  the  Issue  of  target  glint,  we  were  faced  with  the 
question  of  what  to  use  as  a glint  In  the  experimental  configuration, 
and  what  theoretical  model  should  be  used  to  represent  this  glint.  The 
question  arises  because  of  both  the  experimental  constraints  Imposed  by 
the  targets  of  interest  and  by  the  numerical  evaluation  of  the  complex 
theoretical  expressions  that  accurately  model  the  propagation  geometry. 

The  theoretical  model  chosen  is  discussed  in  Section  IIA.  Experimentally, 
we  wanted  to  ensure  that  the  Irradlance  received  was  not  dominated  by 
either  the  diffuse  term  or  the  coherent  (glint)  term  alone,  in  order 
that  the  experimental  results  Included  a realm  where  both  terms  were 
effective  in  influencing  the  statistics  of  the  received  irradlance.  Since 
the  diffuse  portion  of  the  target  scatters  the  incident  radiation  into 
large  angles,  the  received  Irradlance  due  to  the  diffuse  target  is  small. 
This  required  that  the  received  Irradlance  due  to  the  glint  be  small  also, 
and  consequently  that  the  coherent  feature  on  the  target  be  small.  We 
felt  that  a retro-reflector  glint  did  not  satisfactorily  represent  a 
real  glint,  and  since  the  required  radius  of  curvature  of  a curved  specular 
surface  that  would  give  an  appropriate  level  of  irradlance  at  the  target 
was  extremely  large,  the  experimental  glint  was  chosen  to  be  a small  flat 
mirror.  This  was  quite  practical  experimentally  and  gave  results  in  good 
agreement  with  the  theoretical  model  chosen. 

In  addressing  the  issue  of  source  coherence,  we  utilized  an  Argon  ion 
laser  operating  on  the  488  nanometer  line  and  contrasted  the  results  with 
and  without  a temperature  stabilized,  intra-cavity  etalon.  With  the  etalon 
in  place,  the  laser  output  consisted  of  a single  longitudinal  TEM  mode. 


With  the  etalon  removed,  multiple  longitudinal  modes  were  allowed  to  run 
covering  the  full  488  nanometer  llnewldth.  The  llnewldth  Is  on  the  order 
of  10^®  Hertz  and  this  configuration  reduced  the  coherence  length  to 
approximately  3 centimeters. 

The  basic  experimental  configuration  Is  shown  In  Figure  1.  The 
source  Is  an  Argon  Ion  laser  operating  at  488  nanometers  that  can  be 
conveniently  modified  to  run  either  single  or  multimode.  The  propagation 
path  is  nominally  500  meters  over  flat  featureless  farm  land.  The  receiver 
consists  of  two  photomultiplier  tubes,  each  with  a 1 millimeter  aperture. 
The  spacing  of  the  two  apertures  is  variable  to  allow  measurement  of  the 
spatial  correlation  function. 

This  report  Includes  the  results  of  a detailed  theoretical  investi- 
gation of  speckle  propagation  in  atmospheric  turbulence  using  the  extended 
Huygens  Fresnel  principle.  This  investigation  is  set  forth  in  Section  II. 
The  numerical  evaluation  of  particular  theoretical  expressions  is  dis- 
cussed in  Section  III.  The  experimental  program  is  described  in  Section 
IV  and  exprimental  results  are  presented  and  compared  with  the  numerical 
results  where  appropriate.  A discussion  of  the  results  in  Section  V 
concludes  the  report. 


> 
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II.  Theoretical  Description 


A.  Introduction 

Theoretically,  we  are  interested  in  describing  the  statistical 
12  3 

nature  ’ ’ of  the  field  scattered  from  a target  in  the  geometry  shown 

schematically  in  Figure  1.  In  order  to  do  this,  we  use  the  extended 
4 5 

Huygens  Fresnel  ’ principle  to  compute  the  field  incident  on  the  tar- 
get due  to  a focused  gaussian  beam  propagating  over  an  atmospherically 
turbulent  path.  The  field  scattered  from  the  target  is  then  related  to 
the  incident  field  using  the  known  characteristics  of  the  target.  Finally 
this  scattered  field  is  propagated  back  to  the  receiver  using  the  ex- 
tended Huygens  Fresnel  principle  once  again.  This  gives  an  expression 
for  the  field  at  the  receiver  that  may  be  used  to  compute  the  statistical 
averages  of  interest.  In  this  Section  we  compute  the  mean  irradiance  for 
a coherent  source  and  a diffuse  target  with  multiple  glints  in  part  B. 

In  part  C we  evaluate  the  second  order  moments,  namely  variance  and  co- 
variance  of  irradiance.  The  results  of  part  B are  then  extended  to  in- 
clude reduced  source  coherence  in  part  D. 


B.  Mutual  Coherence  Function  and  Mean  Irradiance 

With  r denoting  a two-dimensional  coordinate  vector  in  the  plane 
of  the  transmitter  aperture,  the  source  amplitude  distribution  for  a single 
mode  laser  transmitter  is 

Uo(?)  = exp  ^ - g rA  (1) 

' 0 ' 

I I 

where  r = | r | , 

a = the  radius  where  the  intensity  is  e"^  of  the  intensity  at  r = 0, 

k°  = ^ 

^ \ ' 

X = wavelength,  and 
F = the  focal  length. 

Using  the  extended  Huygens  Fresnel  principle,  the  field  at  the  target  can  be 
expressed 


/*  r I I ^ 


d?  (2) 


where  U^(p)  = the  field  incident  upon  the  target, 

p = a two-dimensional  coordinate  vector  in  the  plane  of  the 
target, 

L = path  length, 

and  ipjCp,  r)  describes  the  effects  of  the  random  medium  (atmospheric  turbu- 
lence) on  the  propagation  of  a spherical  wave  from  the  source  to  the  target. 


For  a target  composed  of  a diffuse  surface  containing  multiple 


glints,  the  field  scattered  from  the  target  can  be  written 


face,  and  U (p)  is  that  portion  of  the  scattered  field  due  to  the  multiple 


glints.  The  field  scattered  from  the  diffuse  portion  of  the  target  can  be 


where  R is  the  reflectivity  of  the  diffuse  surface.  Assuming  that  a cube 


corner  is  a valid  model  for  a small  glint,  the  field  reflected  from  the 


multiple  glints  may  be  written 


the  complex  strength  of  the  m glint 


the  effective  amplitude  reflection  coefficient  of  the  m glint 


the  position  of  the  m glint,  and 


the  number  of  glints 


In  the  following  development,  we  will  assume  that  for  small  glints,  r (p) 


has  a gauss ian  form,  i.e 


where  Ap  = the  "width"  of  the  m glint 


In  terms  of  the  field  scattered  from  the  target,  the  field  at  the 
receiver  may  be  obtained  by  reapplying  the  extended  Huygens  Fresnel  principle. 


Ur(p)  = exp  (ikl)  J exp  [ 


where  U (p)  = the  field  at  the  receiver, 
r 


p = a two-dimensional  coordinate  vector  in  the  plane  of  the 
receiver, 

and  4/2 (P*  P)  describes  the  effect  of  atmospheric  turbulence  on  the  propagation 
of  a spherical  wave  from  the  target  to  the  receiver. 

At  this  point,  we  can  write  an  expression  for  the  mutual  coherence 

function,  r(pi,  P2). 

r(Pi,  P2)  = <Uj.(Pl)  U*(P2)> 


“ // ^’s  (P2)> 


iklpi  - P2I  iklp2  - P2I 


(8) 

X <exp  h;/2(Pl>  Pi)  +4^2*(P2’  P2)J^ 

where  the  * indicates  the  complex  conjugate  and  the  angular  brackets  < > 
Indicate  ensemble  averaging.  The  mutual  coherence  function  of  the  scattered 
field  contained  in  the  Integrand  of  equation  (8)  can  be  rewritten 

<Ug(Pi)  U*(P2)^  ' + Ug(Pl)]  [^d^°2)  + Ug(P2)]  > 


(9) 


- <U^(Pi)  U^(P2)>  + <U^(Pi)  Ug(p2)> 
+ <U^(pi)  U*(P2)>  + <Ug(Pi)  U*(P2)> 


The  diffusely  reflected  fields  are  spatially  incoherent  so  that 


<Uj(Pi)  U^(p2)>  “ '5(pj  - P2) 


where  <Ij(Pi)>  * '^Ujj(Pz)  '^^Pl  “ ^2^  is  the  Dirac  delta 

function.  Also,  due  to  the  spatial  incoherence  of  the  diffusely  reflected 


fields. 


<U^(Pl)  Ug(P2)>  = <Ug(Pi)  U^(P2)^  ■=  0 


<U^(pl)  Ug(P2)>  = <I^(pl)>  <5(p2  - P2) 


+ y\  y]  3 a*  <u*(p  ) u (p  )> 

in^  102  i nij  1 ni2 

mj  = 1 m2  * 1 


|Pl  - P, 


iP2  - Pr 


Following  this  we  will  assume  that  the  glints  are  identical,  i.e.  that 

a “1  and  Ap  * Ap  . 
m^  m^  w 

Since  from  equation  (11),  the  cross  correlation  terms  are  zero, 
the  mutual  coherence  function  of  equation  (8)  can  be  expressed  as  the  sum 
of  a diffuse  component  and  a "specular"  component  (from  the  multiple  glints); 


L 


4- 


r 


Substituting  equation  (1)  Into  equation  (17)  and  using  the  equivalent 
of  equation  (15)  gives 


- o _ 

• " ["  2l  (^  ■ ?) 

* exp  ^ Pz  ■ (Xl  ' X2)j 

r /!"*  * 

/|X1-X2I\  , 


where  “ |r^|.  Substituting  this  expression  Into  equation  (16)  gives  a 
general  result  for  the  diffuse  component  of  the  mutual  coherence  function. 

The  resulting  expression  can  be  evaluated  for  the  focused  case  (F  ■ L)  giving 

r^(Pi»  Pa)  “ 


I Pi  - Pal 


Pi  - Pa 


Using  once  again  equations  (8)  and  (12),  the  specular  component 


(PI.  P2)  - (j^)  ^ Z S /[  '"i(".,)'’l(%2)* 

mi  » 1 m2  “ 1 


•0 


■*  I-*- 

* exp  ^ l^iPl  ■ Pi  I - IP2  ~ P2 1 / 

V2(Pi»  Pi)  + >^2(P2*  P2)J^  ^20 

after  having  incorporated  the  assumptions  given  following  equation  (12). 
Noting  that 

<exp  rJ2(Pl>  Pi)  'l'2(P2>  P2)J^  ' 

exp  j^-  D^i  - P2,  Pi  - P2)] 

where  d(pi  - P2>  Pi  “ P2)  wave  structure  function,  equation  (20) 


can  be  rewritten 


P2)  '(^)  “p[ir(pi'  -P2')] 


M M 


« E E "i(“«2)’ 

X jj  dPidp2  exp  I - -J—  (|P1  - P^J  + |P2  - Pn,J  ) 
+ ^ [(Pi^  - 2pi  • Pi)  - (P2^  - 2P2  • P2)] 


d(p1  - P2»  Pi  - P2).' 


If  we  make  the  change  of  variables 


p “ Pi  “ P2>  snd 

2R  - Pi  + P2* 

the  integration  over  R can  be  performed  giving 


12 


where  p * P]  - P2  and  p « |pj. 

In  the  above  developments,  it  has  been  implicitly  assumed  that 
the  "size"  of  the  glints,  is  much  smaller  than  the  beam  at  the  tar- 
get. If  it  can  be  assumed  further  that  (Ad^)  is  small  compared  with  and 
p^,  then  in  integral  expressions  such  as  Eq.  (20)  all  terms  are  varying  ex- 
tremely slowly  with  respect  to  p compared  to  the  magnitude  of  the  glint  term 
and  the  approximation 


exp 


^ IK  - ir  - D-  '^1 

-^(!r 

■ °(P*  Pm  “% 

\ mi  m^/k 


p - P "Pi 
mj  m2 


(25) 


mi  m2/^ 

Thus,  the  mutual  coherence  function  of  equation  (8)  can  be  evaluated  by  using 
equations  (19)  and  (23)  or  (25)  in  equation  (13). 

These  expressions  can  be  used  to  evaluate  the  mean  irradiance  at 
the  receiver.  For  the  focused  case, 


<I  (p)>  = r(p,  p) 


= rIu  |2  ^ + 


M M 


° l2 


(i)  E E <<(?.,) 

mi  * 1 m2  * 1 


Jpdp  Jo|p|r(p  - “?) 


+ i 


in- 


exp 


where  the  change  of  variables 


( ^ 1.  - (» j ( 


(26) 


* 

0 

- P » 

m- 

mi 

m2 

-¥ 

-► 

D • 

P 

m- 

+ P 

ni2 

has  been  used  and  the  angular  integration  performed.  Using  the  approximation 
of  equation  (24),  this  can  be  written  in  the  simplified  form 

a ^ 

<I(p)>  - R'U  |2 

° L2 
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C.  Covariance  of  Irradlance 

Many  of  the  significant  statistical  properties  of  the  received 
field  are  revealed  by  an  examination  of  the  covariance  of  Irradlance  at 
the  receiver. 

Cj(Pi,  P2)  = <Ij.(Pi)  - <I^(Pi)><Ij.(P2)>  (30) 

In  order  to  evaluate  the  covariance  of  irradlance,  it  is  convenient  to 
focus  on  each  of  the  two  terms  on  the  right  hand  side  of  equation  (30) 
separately.  The  first  term  is  the  correlation  function  of  irradlance: 
Bj(Pi,  P2)  = <I^(P:)  Ij.(P2)^ 

= <U^(Pi)U*(pi)  U^(p2)  U*(p2)>  . (31) 

An  expression  for  this  correlation  function  can  be  written  by  applying  the 
extended  Huygens-Fresnel  principle  to  the  fields  scattered  from  the  target. 
Thus : 


:)//// -=. 


dp2d03dP4  <U^(pi)  U^(p2)  U^(P3)  U^(P4)> 


f 

r 


lie 

2L 


[1;. 


Pi 


IPl  - P2I 


IP2  - P3I 


|P2  - P4i 


1»  Pj)  + P2)  "*■  P3)  ^2(P2>  P4)J^ 

(32) 

Using  equations  (3)  and  (11),  the  fourth  order  correlation  of  the  scattered 
field  in  the  integrand  of  equation  (32)  can  be  rewritten  as 


Ug(oi)  U^(02)  U (03)  U^(04)>  = <U^(pi)  U^(P2)  U^(03)  U^(P4)> 


+ <U  (pi)  U (P2)  U (P3)  U (P4)> 
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+ <Ug(Pi)  U*(P2)  U^(P3)  U*(P4)> 

+ <Ug(pi)  U*(P2)  U^(P3)  U*(P4)> 

+ <U,(Pi)  U*(P2)  U (P3)  U*(P4)> 

+ <U  (Pj)  U*(P2)  U (p,)  U*(P4)> 

g ^ g ^ g ^ g 


^ 'j  A^(pj,  p2»  P3>  P4) 
n = 1 


(33) 


where  the  are  defined  by  equation  (33).  Then  the  correlation  of  irradiance 
at  the  receiver  can  be  written  as  the  sum  of  six  terms: 


P2)  ° P2) 

n = 1 

where  from  equations  (32)  and  (33), 


(34) 


(Pi.  P2)  = fe)  ////  do idp2dp 3dp4  A^(pi,  P2 > P3>  P4) 


exp  i ^ I Pi  - Pi  1 - I Pi  - P2  I 


2L 


+ |P2  - P3I  - |P2  - P4I  ] ^ 


X <exp 


U2(Pl>  P])  + 'i'2  (Pl’  P2)  + ’+'2(P2»  P3)  ’^2 


(P2»  P4)] 
(35) 

The  fourth  order  mutual  coherence  function  due  to  turbulence  in  equation 
(32)  is  defined  as 

H(pi,  P2;  Pl>  P2’  P3*  P14)  ' ^^^P  L^2(Pl’  Pi)  + ^2  (Pl»  P2) 

+ :i'2(P2»  P3)  'i'2(P2’  P4)J^ 

= exp  j - I [d  (Pl  - P2)  - Dg,(Pl  - P2«  Pi  -P3> 
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^ ‘ 


I 


+ D^^Cpi  - P2,  Pi  - P14)  + D,,(Pi  - P2.  Pz  - P3) 


'P 


%(Pl  - P2.  P2  - P4)  + D^(P3  - P4)] 


>1^ 


+ 2Cjj^(pi  - P2,  Pj  - P3)  + 2Cj^(pi  - P2,  P2  - P4) 

(36) 

where  D and  C are  the  wave  structure  function  and  the  log  amplitude  co- 

y'  Jw 

variance  function  for  a spherical  wave,  respectively. 

Using  the  notation  of  equation  (35)  along  with  equations  (10)  and 
(11),  the  six  terms  making  up  the  correlation  function  of  irradiance  at  the 
receiver  may  be  written.  Using 

Ai(Pi.  P2.  P3.  P4)  = (pi)  U (P2)  U (P3)  U (P4)> 


d 

+ <U^(Pi)  U^(p4)><U^(P2)  U^(p3)> 


= <U  (Pi)  U*(?2)><U  (P3)  U*(P4)> 


4it 

.1,2 


<I^(pi )><I^(P3)>  6(Pi  - P2)  5(P3  - P4) 


+ <I^(pi)><I^(P3)>  6(pi  - P4)  6(p3  - P2) 


in  equation  (35)  gives 

B,  (Pi,  P2)  = / 1 ' 

I exp  [^Cj^(pi  - P2>  P?  - P4)]  + exp  (pi  - P2) 


dp2dp4  <I^(p2)><I^(P4)> 


B(PiJ  P2 * P4»  P2>  P2 > P4) 


(02 

(37) 


- P4)j 


-►  -► 


Using  A2(Pi,  P2.  P3.  P4)  = (pi)  U,  (P2)  U (P3)  U (04)> 
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= <U,(pi)  U*(C2)><U  (P3)  U*(P4)> 

d ^ d g g 

= — <I  (Pi)><U  (P3)  U*(p4)>  6(pi  - P2) 

, 2 “ 8 8 


in  equation  (35)  gives 


P2)  = — fjf  dP2dp3dP4  <I^(P2)><Ug(P3)  Ug(P4)> 

X exp  (P3  + P4  - 2p2)  * (P3  - P4)J 

^ “V  . 

^ H(pi,  P2  > P2«  P2» 

Using  A3(Pi,  P2.  P3.  P4)  = ^^^(Pl)  ^’g(P2)  U^(p3)  U^(p4)> 

= <U  (pi)  U (P2)^^^'j(P3)  Ui(P4)> 
g 8 d a 

= — <U  (Pi)  U*(p2)><I,(P3)>  <5(P3  " P4) 
1^2  g ^ 8 d = 


in  equation  (35)  gives 


JJJ"  doidp2dp4  <I^(P4)><Ug(pi)  Ug(P2)> 

[ik  “^  \ \~l 

(Pl  + P2  " 2pi)  • (Pi  - P2)J 
^ H(pi , P2 ; Pi.  P2 » P4>  P4) 


Using  A4(pi,  P2,  P3,  P4)  = 


<Ug(;i)  U*(P2)  U^(P3)  U*(;4)> 

<U  (Pi)  U*(p4)><U*(?2)  U^(P3)> 
S 6 


— <U  (Pi)  U*(p4)><I  ,(p3)>  5(p3  - P2) 
j^2  8 8 d 


in  equation  (35)  gives 
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<I  (p)>  = ^ I (p)> 

r d g 

= <I  (p)>  + <I  (p)> 

d g 

where  <I  (p)>  and  <I  (p)>  may  be  obtained  by  setting  Pi  = P2  “ 
d g 

(19)  and  (22)  respectively.  Using  equation  (43)  then, 

= <I^(Pl)><Ij(P2)>  + 

+ <I  (pi)><I  (P2)>  + (Pi)^*^!  (P2)^ 

g S g g 

An  expression  for  the  covariance  of  irradiance  can  be 
the  following  definitions: 

Cj(Pi,  P2)^  = <Ijj(Pl)  ~ 

= Bj^(Pi>  P2)  “ <Ijj(Pi)><Ij(P2)>> 

Cj(Pi,  P2)g  = ^lg(Pl)  Ig';P2J^  - <Ig<P*l)>^Ig(p2)'’ 

= Bj^(Pi,  P2)  - <Ig(Pl)><Ig(P2)>’ 

Ci^j  (Pl.  P2)  = Ig(P2)>  - 

= Bj^(pi,  P2)  - <Ij(Pi)><Ig(P2)>. 

Cj  j (P2>  P^)  = 

d g ^ ® 

= Bj^(pi,  P2)  - <I^(p2)><Ig(Pi)>. 

Cj,  (Pl.  P2)  = (pi)  U*(pi)  U*(P2)  U (P2)> 


(43) 

in  equations 

(44) 

written  using 

(45) 

(46) 

(47) 

(48) 


= B (pi , P2)  , 


(49) 


and 


<U  (pj)  U*(P2)  U (P3)  U*(?4)> 
8 ^ 8 ^ 8 8 


M M M M 


■ E E E E \{K,)  <(»,3) 

- 1 1112  ~ 1 ®3  ~ 1 ®4  ” 1 


exp 


p.  - " (^3  - ?.3)^ 

* - v)^ 


(50) 


where 


"i(s)  <(“"3) 


dridr2dr3dr4  U (rj)  U(r2)  U (^3)  11(^4) 


■ (*)//// 


jdc 

2L 


in4 


- r. 


(->■  -»•  ->•  ->• 

Pn,  » Pn,  • ’^2*  • 

mi  m2  m3  mu  1 ^ *♦/ 


m^  m2  m3  m4 
The  covariance  of  irradiance  becomes 


(51) 


Ci(Pl.  P2)  = Cj(Pi,  P2)^  + Cj(Pl.  P2)  + C - ) 

d g ^ 


+ C (p2,  P])  + C (Pj,  P2)  + C (Pl.  P2) 
d g ^dg  ^dg 


(52) 


The  variance  of  irradiance  can  be  obtained  from  equation  (52)  by  setting 
Pi  “ P2  ~ P 
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D.  Speckle  Patterns  With  Reduced  Source  Coherence 

8,9  , 

In  this  section,  the  effect  of  reduced  temporal  coherence  ot 
the  source  is  considered.  Following  the  development  of  the  previous  sec- 
tions, and  using  the  definitions  given  there,  the  two  frequency  mutual  in- 
tensity function  is 

rr(Pi>  P2>  ^^2^  <Uj.(Pl>  kj)  U^(P2,  k2)> 


= — ff  dpidp2  < (Pi , kj) 

(2tiL)2  JJ  I ^ 

+ 7 U.(c  k,  ) exp / 

^ A "1.  A Y ^ 

mi  = l L\w/  J; 


-r1/2  U*  (P2.  k2)  + ^ k2) 

m2  = 1 


(11^1/ 


exp  - 


:xp  (Pl  - Pi)  (P2  " Pi) 

:exp|^!|/(p  j , Pl,  ki)  -I-  i|J  (P2  5 P2>  ^2 )] ' 


, ^ _ 2Tt 

where  k,  = -^  , 
i 


= wavelength. 


I 


In  order  to  simplify  the  expression  in  equation  (54),  the 


I 


ensemble  averages  implicit  in  the  integral  can  be  evaluated.  For  the 
diffusely  scattered  radiation. 


4ttR 


<U^(Pl,  ki)  U^(P2,  k2)>  = — 6 ki)>  6(Pi  - P2)  6(ki  - k2)  (55) 


where  6 is  a normalization  constant  that  will  be  evaluated  below.  It  has 


been  suggested  that  the  log  amplitude  covariance  function  and  the  wave 

10 


structure  vary  weakly  with  wave  number.  We  extend  this  to  the  narrow  band 
spherical  wave  for  the  purposes  of  this  calculation  and  assume  that 


<exp[v(pi,  Pj,  ki)  + ip  (P2,  P2.  k2)]> 
can  be  approximated  by 


:exp  j^iii(Pl , pj,  k^)  + li)  (p2,  P2,  k^)J  > = exp[^-  y D(p2  - Pj,  P2  “ Pi,  kj  (56) 


where  k is  the  mean  wave  number.  For  the  coherently  scattered  radiation 

o 


(from  the  glints),  we  have  for  the  focused  case 


<U  ( p , kA  U fp  k2^  > 

8V®1  / SV  "’2«  / 


kik2  ‘ 


r'-Pi-fe) 


(2L) 


iki 


ik2  I 

nT  ■ 2ir  “^m^^  ( 


/ 


X I dr  exp 


- i 


L 4a 


1 /-►  \ 
- 7 D(  p - p , r,  k ) 
2 \ m^  0/ 


(57) 
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-li. 


Is  the  diffuse  contribution  to  the  mutual  intensity, 
1 


o kj 


klPl  - k2P2 


and  the  assumption  has  been  made  that 
2 

1 

expl 


Rr  - 


and  where 


dp  **  kik2 

r^f(Pl.  kj;  P2,  k2)  » — ; — 

(2l)2 


M M 


exp  (kiPi^  - k2P2^)] 


m^  1 1 


Tiki 


ik. 


®*P[2r  (Pm'  - - -2L  (^m'  " 


- 2 °Ip2  - Pi’  Pm  ■ Pm  ’ ^^o) 


(60) 
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is  the  contribution  from  the  specular  portion  of  the  target. 


The  mutual  intensity  at  the  receiver  is  the  integral  of  equation 
(58)  over  the  appropriate  spectral  distribution  of  the  source.  Thus 

TfCPi.  P2)  =j^  dkidka  [^df(Pl«  1^1;  P2.  1^2)  + TgjCpi,  kj;  P2,  k2)] 


X S(ki)  S(k2) 


= P2)  + P2) 


If  the  source  has  a gausslan  spectral  density  with  r.m. s.  deviation  W and 


mean  wave  number  k , then 
o 


S(k)  = exp 

)/2ttW 


(k  - k^)^ 
o 


Using  this  assumption,  the  diffuse  contribution  to  the  mutual  intensity 


at  the  receiver  is 


a 2 u ^ 

pdf^pi’ P2>  = ~;= — ®*p  " 

2/irw 


W2(P1^  - P2^)  |P2  - Pll 


h - Pi  ik 


+ (Pi^  - P2^)  (63) 


Comparing  this  with  previous  results  in  the  limit  of  W ->■  0,  the  nomallzatlon 
constant  can  be  evaluated. 

6 * 2/TrW  (64) 

In  order  to  simplify  the  evaluation  of  the  glint  contribution  to 
the  mutual  intensity,  the  parameters  associated  with  the  Argon  ion  laser 
source,  namely  k^  ■ 1.29  x 10^  , Ak  • 209  and  consequently  ^ ■ 1.63  x 10“^, 
are  considered.  Making  the  change  of  variables,  kj  ■ kj  and  k2  ■ kj  + Ak, 
the  following  approximations  are  useful: 


and 


where 


P„  “ P - P , and 

“21  "*2  “1 


Following  this  same  approach,  expressions  for  the  mean  and  related 
second  moments  can  be  evaluated,  where 

<I(p)>  - J S(k)  <I(p,k)>  dk  (66] 

<I^(P)>  - jfjf  S(ki)  S(k2)  <I(p,ki)I(p,k2)>  dkidk2  (67) 

<I(Pl.P2)>  “jfjf  S(ki)  S(k2)  <I(pi,ki)I(p2,k2)>  dkidk2  (68' 

As  in  the  case  of  the  coherent,  or  monochromatic,  source,  these  expressions 
must  be  evaluated  numerically. 
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III.  Numerical  Evaluation 


In  addition  to  the  numerical  techniques  described  previously, 

Che  following  additional  methods  were  used  to  evaluate  the  expressions 
developed  in  Section  11: 

* The  scotchllte  was  assumed  to  be  equivalent  to  a diffuse 
surface  with  a gain  of  1000. 

• In  the  evaluation  of  ^ in  Eq.  53,  we  used  experimental 

8 12 
data  for  the  variance  of  a focused  gausslan  beam. 

The  principal  results  are  associated  with  the  normalized  variance  of 

irradlance  for  a coherent  source  and  a target  consisting  of  a diffuse 

surface  with  a single  glint.  In  particular,  it  has  been  determined  that 

the  normalized  variance  of  irradlance  can  be  on  the  order  of  3.  This 

corresponds  to  fluctuations  in  the  received  irradlance  that  are  greatly 

enhanced  over  those  present  for  a diffuse  target  alone.  These  enhanced 

fluctuations  may  have  a significant  impact  on  systems  operating  in  this 

geometry. 

Table  1 is  a typical  example  of  the  results  of  these  calculations, 
showing  a peak  in  the  normalized  variance  of  irradlance  of  3.3.  The  param- 
eters used  in  this  calculation  are: 

L “ 500  m (path  length) , 

A “ 488  nm  (wavelength), 

a “1.35  cm  (radius  of  transmitter  beam), 

0 

Og  • 0.19  mm  (radius  of  glint), 

and  both  the  glint  and  the  received  irradlance  are  located  on  the  optical 

axis. 
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TABLE  I 


NORMALIZED  VARIANCE  OF  IRRADIANCE  (a  2) 

^N 

DIFFUSE  TARGET  WITH  SINGLE  GLINT 


<I  > 
g- 


5 X 10” 

,090075 

25.895 

10-15 

.14643 

24.685 

2 X 10-15 

.35566 

22,542 

5 X 10”15 

1.0767 

17.742 

10-1** 

1.9127 

12.909 

1.5  X 10”  1** 

2.  6526 

10.059 

2 X 10”  1** 

3,3154 

8.434 

3 X lo”!** 

3.1741 

6.0089 

5 X lO"!** 

2. 7002 

3.9666 

10-13 

2.1533 

2.1015 

2 X 10-13 

1.6054 

.72014 

5 X 10-13 

1.1498 

.029616 

10”  12 

1.0921 

. 000145 

refraction 

structure  constant 

<1  > * Mean  Irradiance  from  glint:  <I,>  = Mean  irradiance  from  diffuse 
g d 

portion  of  target 
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Note  that  at  low  turbulence  levels  (small  C„^),  the  glint  dominates 

<I  > ^ 

the  received  irradiance,  i.e.  • ^ >>  1,  while  at  high  turbulence  levels, 

"^d^  <I  > 

the  diffusely  scattered  light  dominates  the  received  irradiance,  <<  1. 

Thus  we  expect  to  see  the  normalized  variance  of  irradiance  return  to  unity 

at  high  integrated  path  turbulence  as  for  a diffuse  target  alone. 


Figures  2 through  6 Illustrate  the  results  of  the  numerical 
evaluation.  Figure  2 shows  the  results  of  previous  work  for  a coherent 
source  and  a diffuse  target  for  purposes  of  reference.  The  solid  curve 
Is  a numerical  result,  and  experimental  measurements  are  superimposed  on 
the  graph.  Several  features  of  this  curve  are  of  Interest.  At  low  tur- 
bulence levels,  0^*1  due  to  the  diffuse  nature  of  the  target.  At 
high  turbulence  levels,  o ^ ■ 1 also.  This  is  the  familiar  saturation 
region,  and  for  this  particular  target  the  normalized  variance  returns 
to  unity.  At  intermediate  levels  of  turbulence,  o ^ peaks  up  at  approxi- 
mately  1.25.  This  peak  is  due  to  the  log  amplitude  terms  in  the  mutual 
coherence  function  as  described  previously. 

Figure  3 Illustrates  the  effect  on  o ^ with  the  addition  of  a 
glint  or  specular  component  to  the  target.  For  this  case,  the  glint  is 
located  1 cm  off  axis;  that  is,  the  glint  is  displaced  by  1 cm  from  the 
center  of  the  focused  gausslan  beam  on  the  target.  (For  no  turbulence 

the  — beam  radius  at  the  target  is  calculated  to  be  0.66  cm).  For  the 

e 

smallest  glint  shown  (0.1  mm  radius),  the  irradiance  at  the  receiver  is 
always  dominated  by  the  diffuse  portion  of  the  target.  Consequently  the 
curve  labelled  0.1  mm  is  practically  identical  to  that  for  a diffuse  tar- 
get alone.  (Note  the  change  of  scale  on  the  ordinate  between  Figures  2 
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Figure  2.  Normalized  variance  of  Irradlance  vs  log  amplitude  variance  for  a 

diffuse  target  Illuminated  by  a coherent  source  over  a 500  meter  path. 
Solid  line  Is  the  theoretical  curve  and  ® represent  experimental  data. 


with  a single  glint  illuinlnated  by  a coherent  source  over  a 500  m path.  Each 
curve  corresponds  to  a different  size  glint  with  the  1/e  radius  of  the  model 
gausslan  glint  given.  The  glint  is  located  1 cm  off  the  optical  axis. 


nt  Illuminated  by  a coherent  source  over  a 500  m path.  Each  curve 
different  size  glint  with  the  1/3  radius  of  the  model  gaussian 
I glint  is  located  0.5  cm  off  the  optical  axis. 


and  3.)  As  the  size  of  the  glint  increases,  however,  the  glint  term  be- 
comes important.  This  causes  two.  effects:  1)  at  low  turbulence  levels. 


n ^ is  reduced  from  unity  as  would  be  expected,  and  2)  a ^ peaks  up  to 
3.8  at  moderate  turbulence  levels.  These  effects  are  evident  in  Figures 
A and  5 also,  where  the  only  change  from  Figure  3 is  that  the  glint  is 
located  0.5  cm  off  axis,  and  on-axis  respectively. 

Figure  6 shows  the  normalized  variance  of  irradiance  vs  a ^ for 

X 

a specular  feature  of  0.19  mm  radius,  for  three  glint  locations;  on  axis, 

0.5  cm  off  axis  and  1 cm  off  axis.  This  curve  will  be  used  for  comparison 
with  experimental  results  in  Section  IV. 

One  final  feature  of  interest  is  evident  in  Figures  3-5.  That 
is,  the  tendency  of  o ^ for  the  largest  glint  (0.5  mm  radius)  to  "hang  up" 
at  a value  near  3 before  falling  to  unity  at  high  turbulence  levels.  This 
will  be  discussed  in  the  Append' x where  a simplified  model  is  suggested 
for  the  large  normalized  variance. 

Numerical  evaluations  of  the  analytical  expressions  derived  in  the 
multiple  glint  case  and  the  partially  coherent  source  case  have  not  been 
completed.  With  current  techniques,  the  multiple  glint  case  takes  im- 
practically  long  to  evaluate  and  is  not  expected  to  contribute  anything 
f undame" tally  different  than  the  simple  glint  case.  With  regard  to  the 
partially  coherent  source,  there  seems  to  be  something  missing  in  the 
analytical  understanding  of  the  diffuse  scattering  that  leads  to  substantial 
differences  between  theory  and  experiment.  This  case  is,  consequently,  not 
well  understood. 


r 


38 


with  a single  glint  illuminated  by  a coherent  source  over  a 500  m path.  The 
curves  are  for  a model  gaussian  glint  with  a 1/e  radius  of  0.19  mm  located 
on  axis,  0.5  cm  off  axis,  and  1 cm  off  axis. 


" — 

I 

F- 

IV.  Experimental  Results 

i. 

A.  Introduction  '■ 

An  exprimental  program  to  verify  the  theory  presented  in  Section 
II  has  been  conducted.  In  concept,  the  experiments  were  designed  to  in- 
vestigate two  separate  effects;  one  due  to  target  glint  and  one  due  to 
reduced  source  coherence.  The  experimental  measurements  were  made  using 
the  facilities  and  instrumentation  described  previously. 

A Coherent  Radiation  Model  52  Argon  ion  laser  operating  at  A88 
nm  is  used  as  the  source.  The  laser  has  an  output  power  on  the  order  of 
1 watt  and  contains  an  intra-cavity  etalon  that  can  be  used  to  select  a 
single  longitudinal  mode  when  desired.  In  the  transmitter,  the  laser 
light  is  amplitude  modulated  at  100  k.Hz  using  a Morvue  PEM-e  photo-elastic 
modulator  and  a calcite  polarizer.  The  off-axis  beam  from  the  polarizer 
is  blocked  and  the  scattered  light  is  detected  by  a silicon  photodiode  to 
provide  a 100  kHz  local  oscillator  for  use  in  the  synchronous  demodulators. 

The  on-axis  beam  is  passed  through  a convex  focussing  lens,  a spatial 
filter,  and  a 7 . 5 cm  diameter  objective  lens.  The  resultant  beam  is 
steered  to  the  target  by  a large  flat  mirror.  The  1/e  point  of  the  in- 
tensity profile  is  1.35  cm  at  the  output  mirror. 

The  receiver  contains  two  detectors  mounted  such  that  measurements 
can  be  made  at  various  detector  separations.  Each  detector  consists 

t ^ 

of  an  S-20  photomultiplier  tube  isolated  from  excess  background  light 
by  a 10  nm  optical  bandpass  filter,  and  field  of  view  limiting  apertures. 

The  output  of  the  photomultiplier  tube  feeds  an  operational  amplifier  ' 

and  the  resulting  signal  is  passed  through  a bandpass  filter  to  the 
synchronous  demodulator.  The  demodulator  output  is  low-pass  filtered 
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and  then  recorded  on  a Hewlett-Packard  Model  HP-3960  fm  tape  recorder 


for  subsequent  processing. 

B.  Target  Glint  j 

( 

The  target  used  in  this  portion  of  the  experiment  was  designed  1 

to  provide  a specular  feature  or  "glint"  on  a diffuse  reflecting  surface.  | 

As  in  previous  work  the  diffuse  surface  consisted  of  Scotchlite  (3M  j 

Sprint  Marking  paper).  Substantial  thought  was  given  to  the  implemen-  ! 

tation  of  the  "glint"  portion  of  the  target.  The  eventual  decision  to  ' 

use  a flat  mirror  for  the  majority  of  the  experimental  work  was  based  I 

upon  the  following  reasoning: 

A "real"  glint  consists  of  a specular  reflection  from  a convex  j 

curved  surface.  Experimentally  then,  our  glint  should  consist  of  such  a | 

surface  as  opposed  to  e.g.,  a retroref lector . Additionally,  for  this 
experiment,  we  wanted  the  glint  to  be  such  that  the  field  strength 
at  the  receiver  contained  comparable  contributions  from  both  the  glint 
and  from  the  diffuse  surface.  Otherwise  one  or  the  other  component  would 
dominate  the  statistics.  In  order  to  satisfy  these  two  conditions,  the 
required  radius  of  curvature  for  the  convex  specular  surface  glint  was 

so  large  that  a flat  mirror  provided  the  practical  equivalent.  ' 

t' 

In  order  to  compare  the  experimental  data  with  the  analytical 
predictions,  we  needed  to  compare  our  analytical  model  of  a glint  with  | 

the  flat  mirror  actually  used.  The  analytical  model  was  essentially  that  i ^ 

i '! 

of  a flat  mirror  with  a gausslan  reflectance  distribution.  In  order  to 

compare  this  with  the  experimental  glint  (a  flat  mirror  with  a uniform  ' i 

ii 

I 

reflectance  distribution)  the  1/e  point  of  the  reflectance  was  chosen  to 
give  the  same  on-axls  mean  intensity  as  that  expected  from  the  actual 

I 
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glint  used.  Using  this  as  a basis,  we  obtained  good  agreement  between 
analysis  and  experiment. 

The  variance  results  are  illustrated  in  Figures  7 and  8.  Figure 
7 shows  the  normalized  variance  of  irradiance  vs  log  amplitude  variance 
for  a 0.5  mm  diameter  glint  (the  analytical  model  is  a flat  glint  with 
a gaussian  reflectance  distribution  with  a 1/e  radius  of  0.19  mm).  Experi- 
mental points  are  shown  as  well  as  the  three  analytical  curves  of  Figure  6. 
The  curves  correspond  to  a glint  located  on  axis,  0.5  cm  off  axis  and  1 
cm  off  axis.  These  three  curves  are  shown  to  illustrate  the  potential 
variability  of  the  data  with  respect  to  pointing  the  transmitted  beam 
at  the  glint  and  maintaining  the  pointing  throughout  the  experiment.  Note 


in  particular  that  normalized  variances  as  high  as  2.4  were  measured  in  this 
experiment . 

Figure  8 shows  experimental  measurements  and  analytical  curves 
for  a 4.76  mm  diameter,  537  mm  radius  of  curvature  spherical  glint 
(modeled  by  a 0.144  mm  1/e  radius  gaussian  glint).  Note  again  the  agree- 
ment and  in  particular  the  measurement  of  normalized  variances  on  the 


order  of  1.5. 

In  addition  to  variance  measurement,  the  receiver  was  designed 
to  measure  the  covariance  of  irradiance.  Selected  covariance  curves  are 


shown  in  Figures  9-12. 

Figure  9 shows  the  behavior  of  the  covariance  of  Irradiance  as 
a function  of  the  "strength"  of  the  glint,  with  the  turbulence  level 
constant.  (By  glint  strength  we  mean  the  magnitude  of  the  mean  irradiance 
at  the  receiver  which,  for  a flat  glint,  is  directly  proportional  to  the 


size  of  the  glint  at  the  target.)  As  can  be  seen  in  the  figure,  for  a 
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Normal i zed  variance  of  irradlance  vs  log  amplitude  variance  for  a diffuse  target 
with  a single  glint  illuminated  by  a coherent  source  over  a 500  m path.  The 
curves  are  for  a model  gausslan  glint  with  a 1/e  radius  of  0.19  mm  located  on 
axis,  0.5  cm  off  axis,  and  1 cm  off  axis.  Experimental  points  are  indicated 
on  the  figure. 


Normalized  variance  of  Irradlance  vs  log  amplitude  variance  for  a diffuse 
target  with  a single  glint  Illuminated  by  a coherent  source  over  a 500  m 
path.  The  curves  are  for  a model  gausslan  glint  with  a 1/e  radius  of 
0.144  mn  located  on  axis,  0.5  cm  off  axis,  and  1 cm  off  axis.  Experimental 
points  are  Indicated  on  the  figure. 


Figure  9.  Covariance  of  irradiance  vs  separation  (d)  as  a function 
of  glint  strength  for  a fixed  turbulence  level.  Solid 
line  - no  glint  (diffuse  target  only),  solid  circles  - 
spherical  glint,  open  squares  - 0.48  mm  diameter  flat 
glint,  open  triangles  - 1.00  mm  diameter  flat  glint. 

Path  length  ■ 500  meters,  o ^ “ 0.85. 


fixed  separation,  the  covariance  increases  with  glint  strength. 

Figure  10  shows,  for  a fixed  glint,  the  narrowing  of  the  co- 
variance  with  increasing  turbulence..  Note,  however,  the  "tails"  on  the 
covariance  curves.  These  are  examined  more  carefully  in  the  next  two 
figures,  where  the  effect  is  more  evident. 

To  better  understand  the  shape  of  the  covariance  curves,  con- 
sider the  effect  of  beam  wander  in  the  focused  gaussian  beam  illuminating 
12 

the  target.  As  the  Incident  beam  moves  around  on  the  target  due  to 
beam  wander,  the  illumination  of  the  glint  (which  in  our  case  is  small 
compared  to  the  size  of  the  beam)  changes  and  from  the  point  of  view  of 
the  receiver,  the  glint  appears  to  flash  off  and  on.  Thus  if  beam  wander 
were  the  only  source  of  intensity  fluctuations,  the  irradiance  distri- 
bution at  the  receiver  would  change  uniformly  with  the  illumination  of 
the  glint  and  the  covariance  curve  would  be  flat  with  a value  of  unity, 
independent  of  the  separation.  Under  conditions  where  beam  wander  is 
important  and  where  the  glint  appears  bright  relative  to  the  diffuse 
portion  of  the  target,  we  might  expect  a broad  covariance  curve.  This 
is  shown  in  the  top  curve  of  Figure  11  where  the  glint  is  a flat  specu- 
lar surface  of  2.38  mir  diameter  and  where  D/p  is  near  unity.  Conse- 

0 

quently,  with  D/p  = 1,  beam  wander  is  important  and  we  see  a broad  co- 

o 

variance  curve.  The  lower  curve  is  for  D/p^  = 0.2  with  approximately 
the  same  log  amplitude  variance  0.025). 

This  same  effect  is  illustrated  in  Figure  12.  Under  conditions 
where  the  effect  of  both  wander  and  scintillations  are  evident  in  the 
covariance  curve,  the  covariance  will  appear  to  saturate  at  a value  that 
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is  proportional  to  the  strength  of  the  glint.  This  is  evident  in  Figure 
12  where  the  covariance  saturates  at  a higher  value  for  the  strongest 


glint.  These  covariance  curves  have  not  been  compared  with  analytical 
results  due  to  the  complexity  of  the  numerical  evaluation  of  the  co- 
variance  expressions. 

C.  Partially  Coherent  Source 

In  order  to  investigate  the  effect  of  temporal  coherence  on 
the  statistics  of  speckle  propagation  in  turbulence,  we  repeated  the 
experiments  using  a completely  diffuse  target  (no  glint)  composed  of 
scotchlite,  and  allowed  multiple  longitudinal  modes  to  oscillate  in 
the  argon-ion  laser  source.  This  was  then  compared  with  previous 
results  obtained  using  a completely  diffuse  target  and  a coherent 
source  (single  longitudinal  mode)  as  shown,  e.g.,  in  Figure  2.  We 

; 

found  unexpectedly  that  for  the  partially  coherent  source,  the  normalized  . ■ 

variance  of  irradiance  was  reduced  substantially  from  unity  under  all 

' i 

conditions.  A typical  example  is  shown  in  Figure  13. 

This  result  was  unexpected  because  the  bandwidth  of  the  multi- 
mode  source  was  less  than  or  equal  to  10^*^  Hz,  corresponding  to  a co-  ji 

herence  length  of  approximately  3 cm.  Analytically,  we  expected  no  -jl 

change  (less  than  one  percent)  in  the  normalized  variance  of  irradiance 

' 'i 

for  these  conditions.  The  same  confusion  exists  in  the  absence  of  tur-  j *1 

i i 

bulence.  That  is,  for  zero  turbulence,  we  expect  that  for  the  band-  j i! 

width  given  above,  the  normalized  variance  of  irradiance  will  be  0.999  * ^ 

for  a diffuse  surface.  When  we  perform  an  experiment  in  the  laboratory, 

I 

over  short  path  lengths  such  that  turbulence  may  be  neglected,  we  j 

measured  normalized  variances  on  the  order  of  0.3.  This  was  found  to  ( 
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Figure  13.  Normalized  variance  of  Irradlance  vs  log  amplitude  standard 
deviation  for  a diffuse  target  illuminated  by  a partially 
coherent  source. 


be  true  with  scotchlite  targets,  and  also  for  a variety  of  other  tar- 
gets simulating  a diffuse  surface. 


Since  we  were  unable  to  understand  the  mechanism  for  the 
change  in  speckle  statistics  (which  appeared  unrelated  to  turbulence 
effects)  within  the  scope  of  this  program,  we  concentrated  our  efforts 
on  the  target  glint  measurements  already  discussed. 

i 


( 

I 


V. 


Discussion 


This  program  has  led  to  several  interesting  results  that  include: 

• The  normalized  variance  of  the  irradiance  scattered  from 
a target  containing  a specular  feature,  and  illuminated 
by  a coherent  source  over  a turbulent  atmospheric  path, 
may  reach  a value  on  the  order  of  3 (see  Appendix  for 
physical  model) . This  corresponds  to  fluctuations  that 
are  much  larger  than  for  a diffuse  target  (normalized 
variance  on  the  order  of  unity)  and  will  have  a sub- 
stantial Impact  on  optical  systems  operating  in  this 
geometry. 

• The  normalized  variance  of  the  irradiance  scattered  from 
a diffuse  target  illuminated  by  a source  with  reduced 
temporal  coherence  is  substantially  less  than  when  a co- 
herent source  is  used,  and  in  addition,  substantially 
less  than  would  be  expected  theoretically.  This  experi- 
mental result  appears  to  be  independent  of  atmospheric 
effects  and  is  primarily  a source  target  interaction. 

This  result  is  not  well  understood  and  bears  further 
investigation . 

Both  of  these  results  are  significant  in  the  sense  that  optical 
systems  operating  in  the  presence  of  atmospheric  turbulence,  e.g.  COAT 
systems,  coherent  radar  systems,  etc.,  must  operate  in  the  presence  of 
these  effects.  A continuation  of  the  investigation  of  speckle  propa- 
gation in  the  presence  of  turbulence  will  lead  to  improved  system 
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design  and  analysis,  and  to  improved  system  performance.  Finally  a 
knowledge  of  the  speckle  propagation  statistics  is  necessary  not  only 
for  the  design  of  active  optical  systems,  but  for  use  in  understanding 
countermeasures  against  these  systems  . 
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A simple  physical  picture  of  the  enhanced  normalized  variance  from 
a single  glint  illuminated  by  a coherent  source  over  an  atmospheric  path 
is  suggested  in  this  appendix.  Consider  a bright  glint  so  that  the  dif- 
fuse portion  of  the  target  can  be  neglected.  Assume  that  the  effects  of 
turbulence  on  the  illuminating  beam  and  on  the  reflected  beam  are  inde- 
pendent. Then,  from  the  point  of  view  of  the  receiver,  the  glint  may 
be  considered  as  a source  whose  magnitude  is  modulated  by  the  scintil- 
lations and  angle-of-arr ival  fluctuations  of  the  illuminating  beam.  Then 
the  magnitude  of  the  source  may  be  modeled  by  a random  variable,  X^,  that 
Is  directly  related  to  the  irradiance  of  the  illuminating  beam  at  the  tar- 
get. The  effect  of  the  turbulence  on  the  return  path  may  be  modeled  as  a 

random  variable,  X . Then  the  random  variable  X that  represents  the  ir- 
K 

radiance  at  the  receiver,  is  the  product  of  Xg  and  X^^,  and  the  normalized 
variance  of  X is: 


<x2>  - <X>2  <x-> 

0 ^ 

^ <X>2  <X>2 


1 


1 


I'll  . 1 


Now  consider  the  saturated  case  where  the  normalized  variance  of 
(and  Xg)  is  unity,  i.e. 
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